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ABSTRACT 

We present a detailed X-ray spectral analysis of the Circinus Galaxy using archival 
data obtained with the ASCA satellite. The spectrum shows numerous emission lines 
in the soft X-ray band from highly ionized ions, as well as Compton reflection and 
fluorescent lines from neutral or near-neutral matter. We analyze the spectrum in 
the context of a self-consistent recombination cascade model and find that a nearly 
flat differential emission measure (DEM) distribution in ionization parameter fits the 
data. For a fixed solid angle distribution of matter surrounding a point source, this 
corresponds to a run of electron density of the form, n(r) ~ r- -3 / 2 ^ which is suggestive of 
Bondi accretion onto a central compact mass. Using this density profile and comparing 
the resulting emission spectra with the ASCA data, the size of the X-ray emission line 
region is estimated to be < 1 kpc. We also show that the derived density as a function 
of radius is compatible with the X-ray recombination line emission region being the 
confining medium of the coronal line regions. 

Subject headings: galaxies: individual (Circinus Galaxy) — galaxies: active — galaxies: 
Seyfert — galaxies: starburst — X-rays: galaxies 



1. Introduction 

X-ray spectra of Seyfert 2 galaxies show numerous discrete spectral features that are produced 
in a wide variety of physical conditions. In the few cases where the statistical quality of the data 
has permitted a detailed investigation of the spectrum, the soft X-ray band has been found to 
be almost completely dominated by line emission while the hard X-ray emission has been found 
to be consistent with nonthermal radiation possibly produced through accretion onto the central 
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compact object. X-ray observations have placed meaningful constraints on the structure of matter 
immediately surrouding the central source, and have played an important role in the development 
of a unified model of AGN (Antonucci 1993). Observations of warm absorbers in Seyfert 1 galaxies, 
for example, and their response to variations in the X-ray continuum provide compelling evidence 
that a significant amount of circumnuclear material is photoionized by the primary radiation (Otani 
et al. 1996). There is also observational evidence that a large fraction of Seyfert 2 galaxies har- 
bor starburst regions where large amounts of mechanically-heated gas are produced, for example, 
through supernova explosions. Nuclear jets may also heat the surrounding gas through shocks. The 
orgin of the observed soft X-ray emission and the excitation mechanisms responsible for producing 
the observed emission lines in these galaxies, however, are not well understood and, in most cases, 
cannot be discerned from the low spectral resolution data available to date. 

In the archetypal Seyfert 2 galaxy NGC 1068, for example, Ueno et al. (1994) find that the 
soft X-ray spectrum obtained with ASCA is consistent with emission from two-temperature gas 
in collisional ionization equilibrium (CIE), which may be associated with nuclear and spatially 
extended starburst regions seen in the ROSAT HRI data (Wilson et al. 1992). On the other hand, 
Netzer & Turner (1997) argue that a combination of warm (T ~ 1.5 x 10 5 K) and hot (T ~ 3 x 10 6 K) 
gas in photoionization equilibrium (PIE) is also consistent with the ASCA data. Despite the fact that 
both models adequately fit the data, their are problems in the inferred parameters that cast doubt 
on these interpretations. The CIE fit requires a low-temperature component which is essentially 
metal- free (abundances ~ 3% of the solar photosphere) . The PIE model requires a steep soft X-ray 
continuum component (photon index T = 3.4) that contains a large fraction of the total X-ray 
luminosity and does not have an obvious physical interpretation. On the other hand, a detailed 
X-ray spectral analysis of another Seyfert 2 galaxy, Mkn 3, by Griffiths et al. (1998) shows that 
emission from plasmas both in CIE and PIE can fit the composite Ginga, ASCA, and ROSAT data. 

The Circinus Galaxy is the closest (D ~ 4 Mpc; Freeman et al. 1977) and one of the brightest 
Seyfert 2 galaxies in the X-ray band. Extensive studies at various wavelengths have provided firm 
evidence for an obscured Seyfert nucleus at the center of the galaxy; the detection of a prominent 
one-sided [O III] ionization cone (Marconi et al. 1994), optical and near-IR coronal line emission 4 
from the nuclear region (Oliva et al. 1994; Moorwood et al. 1996), bipolar radio lobes orthogonal 
to the plane of the galaxy (Elmouttie et al. 1995), rapid variations of H2O maser emission from 
the nucleus (Greenhill et al. 1997), polarized and broad Ha emission from the nucleus (Oliva et 
al. 1998), and a highly absorbed hard X-ray spectrum with a bright iron K fluorescent line near 
6.4 keV (Matt et al. 1996; Guainazzi et al. 1999). The system is also known to contain regions of 
starburst activity around the nucleus that account for a large fraction of the bolometric luminosity 
(Rowan- Robinson & Crawford 1989). Recent studies of the infrared spectrum indicate that the 



4 The term coronal lines here refers to high excitation forbidden lines, which were first observed in the solar corona. 
These lines can, in principle, be produced either through photoionization or collisional ionization, and are unrelated to 
emission lines from coronal plasmas, a term which refers to optically thin emission from collisionally ionized plasmas. 
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relative contributions from starburst and AGN activity, which can be traced through emission from 
hot dust and ionized gas, respectively, are approximately equal (Maiolino et al. 1998; Genzel et al. 
1998). 

In a preliminary analysis of the ASCA data by Matt et al. (1996), it was shown that the 2-10 
keV spectrum is dominated by reflection from cold matter, possibly originating from the inner wall 
of the putative molecular torus that surrounds the central nucleus. Several emission lines were also 
detected in the soft X-ray region (0.5 - 2 keV), most of which were inferred to be present in the 
lower resolution BeppoSAX spectrum as well (Guainazzi et al. 1999). The measured energies of these 
emission lines indicate that they are produced through excitation of highly ionized material, and 
cannot be produced in a cold reflection medium. The soft X-ray region is certainly not dominated 
by cold reflection. 

The purpose of this paper is to perform a detailed analysis of the soft X-ray spectrum of the 
Circinus Galaxy using various spectral models, and to derive physical parameters that characterize 
the structure of the X-ray emission line regions. As we show, the soft X-ray spectrum cannot 
distinguish spectroscopically between collisionally ionized gas and photoionized gas (cf., Netzer, 
Turner, & George 1998). There is a strong piece of evidence, however, against an interpretation 
in terms of collisionally ionized gas, and the inferred parameters somewhat favor photoionization 
equilibrium. We discuss the implications of our results and investigate the properties of the active 
nucleus and its circumnuclear environment. 

2. Observation and Data Reduction 

The Circinus Galaxy was observed with ASCA during August 14 - 15, 1995. For the purposes 
of our analysis, we only use data obtained by the Solid-state Imaging Spectrometers (SIS), which 
have the highest spectral resolving power in the soft X-ray band. The SIS data were taken in FAINT 
and BRIGHT 2-CCD modes, while most of the photons from the Circinus Galaxy were recorded on 
SIS0-CHIP1 and SIS1-CHIP3. A quick analysis of the SISO and SIS1 data indicated a ~ 1% offset in 
the energy scale in the iron-K region, while no obvious offset was observed in the soft X-ray band 
(although this is most likely due to the poor statistical quality of data in that region). We applied 
a ~ 1% gain correction to the SIS1-CHIP3 response matrix. Since the spectral resolving power of 
the GIS is lower by a factor of ~ 3 in the soft X-ray region, the discrete features observed in the 
SIS are smeared out in the GIS. We, therefore, only use data collected on the SIS for our spectral 
analysis. 

The data were screened though the standard criteria using FTOOLS v4.0. After conversion 
of FAINT mode data to BRIGHT mode data, we combined them with true BRIGHT mode data to 
obtain a total exposure time of 33 ksec amounting to ~ 4000 and ~ 3000 counts on the SISO and SIS1 
detectors, respectively. Although this method does not fully exploit the resolving power capabilities 
of the SIS, we can significantly improve the statistical quality of the data by compromising for a 20 
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- 30% degredation of the resolving power. 

3. Spectral Analysis 

In this section, we discuss the spectral models that we use for our analysis of the SISO data. 
All of the models we consider here consist of (1) a reflection component, (2) a soft continuum 
component, and (3) a soft X-ray emission line component. Possible physical interpretations of each 
of the spectral components are as follows. The reflection component represents a fraction of the 
intrinsic AGN continuum radiation that is Compton-scattered in a cold medium into our line of 
sight. The spectrum resembles that of a highly absorbed continuum with a number of K-shell 
fluorescent lines superimposed (Lightman & White 1988; George & Fabian 1991; Matt, Brandt, 
& Fabian 1996), which we simply represent by a set of gaussian lines and a Compton reflected 
powerlaw continuum (i.e., PEXRAV in XSPEC). This component has a number of free parameters; 
line energies and fluxes of K-shell Ne, Mg, Si, S, Ar, Ca, and Fe, the intrinsic powerlaw photon index, 
and the total reflected flux. We assume that the K-shell fluorescent lines are infinitely narrow, which 
is justified by the observed Fe Ka line profile, as we discuss later. The soft continuum component 
may be produced through electron scattering of the primary radiation in the surrounding gas, in 
which case the photon index should be equal to that of the primary radiation, or through inverse 
Compton scattering of infrared photons with relativistic electrons. For the latter case, the photon 
index should match the photon index in the radio band. The soft X-ray emission line component 
is modeled as either a collisionally ionized or a photoionized plasma in equilibrium. We assume 
that the entire spectrum is absorbed by a single column density (Galactic + local absorption). The 
reflection component possibly suffers higher absorption, since the reflection region (the inner wall 
of the molecular torus) is assumed to be located closer to the nucleus. The reflected spectrum, 
however, is not sensitive to Njj < 10 22 cm" 2 since most of the flux is contained in the region above 
~ 2 keV anyway. Throughout this paper, we adopt a distance of 4 Mpc to the Circinus Galaxy. 

3.1. Collisional Ionization Equilibrium Model 

We first adopt a model in which the soft X-ray component is represented by a single tem- 
perature thermal plasma in collisional ionization equilibrium (i.e., MEKAL in XSPEC), with the 
temperature, overall abundance, and the emission measure free parameters in the fit. For the re- 
flection component, we first allowed the fluorescent line energies to be free, and found that the 
model is unable to locate the centroid energies for all but the Fe Ka and K/3 lines. We, therefore, 
froze the energies of all of the fluorescent lines, except for those of iron, at values shown in Table 1, 
which are those derived from the photoionization equilibrium fit discussed in § 3.2. These are more 
or less consistent with the ionization state inferred from the iron Ka line energy (i.e., partially 
ionized medium), which is also consistent with the iron Ka line energy measured with BeppoSAX 
(Guainazzi et al. 1999). The measured fluorescent line intensities are also listed in Table 1. Con- 
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trary to Matt et al. (1996), we find a statistically acceptable fit with Xr = 1-21 for 305 degrees 
of freedom. The best-fit model in comparison with the data is shown in Figure 1 and the derived 
parameters are listed in Table 2. 

The allowed ranges in temperature and emission measure (EM = J r? e dV) are shown in 
Figure 2 where we plot 68, 90, and 99% confidence ranges for two interesting parameters. We 
obtain 90% confidence ranges of kT = 0.69±g;gf keV and EM = (1.12+°;^) x 10 63 cm" 3 . Including 
an additional thermal component does not improve the fit significantly (A% 2 ~ 1). In fact, if 
we adopt a powerlaw differential emission measure (DEM) model where the emission measure as 
a function of kT is proportional to (kT) a , we obtain a steep slope of a ~ 5 and a sharp high- 
temperature cutoff at kT ~ 0.7 keV, which nearly corresponds to a single temperature model at 
kT ~ 0.7 keV. Similar results are obtained for different DEM models in coronal equilibrium. In any 
case, the total emission measure in the thermal component is ~ 1 x 10 63 cm -3 and the corresponding 
unabsorbed luminosity is L x = 3 x 10 40 erg s _1 . 



3.2. Photoionization Equilibrium Model 

We next test a model in which the soft X-ray emission lines are produced through cascades 
following radiative recombination in photoionization equilibrium. This spectral model was originally 
constructed for the analysis of the ASCA spectrum of Cygnus X-3 (Liedahl & Paerels 1996) and 
was also succesfully applied to the ASCA data of Vela X-l (Sako et al. 1999). The model includes 
line emission from hydrogen- and helium-like carbon, nitrogen, oxygen, neon, magnesium, silicon, 
sulfur, argon, and calcium, as well as the highest ten charge states of iron (Fe XVII - Fe XXVI), and 
the associated radiative recombination continua (RRC). For a detailed description of the atomic 
model, see Sako et al. (1999). 

We first start by fitting the data with the recombination cascade model, using only the K- 
shell ions of oxygen, neon, magnesium, silicon, sulfur, argon, calcium, and iron, to construct an 
empirical DEM. Due to the limited statistical quality of the data, we cannot directly constrain the 
temperature of each ion emission region from the shape of the RRC, so we fix the temperature of 
each ion at their respective temperature of formation. We do this by first calculating the ionization 
balance and temperature as functions of the ionization parameter (£ = L x /(n p r 2 ), where L x is 
the X-ray luminosity between 1 and 100 Rydbergs, n p is the proton number density, and r is 
the distance from the ionizing source; Tarter, Tucker, & Salpeter 1969) using the photoionization 
code, XSTAR (Kallman & Krolik 1995). We assume that the medium is optically thin with respect 
to continuum absorption, which is an excellent approximation, as we show later. We adopt an 
exponential cutoff powerlaw ionizing continuum with a photon index of T = 1.7 with an e-folding 
energy of Ef = 100 keV and a cutoff energy of E c = 60 keV (Guainazzi et al. 1999). Using 
these results, we then calculate the temperature at which each of the ion line emissivities is a 
maximum and fix each ionic temperature at those values, which we refer to as the temperatures of 
formation. This is a good first order approximation since a large fraction of the total ion line fluxes 
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are produced near those temperatures. Approximately 50% of the recombination line emissivity of 
the He-like ions originates from a range in Alog£ of ~ 0.5. For the H-like ions, the width is broader 
(Alog£ ~ 1) since the emission lines are formed through recombination onto a bare nucleus, which 
exists in a wider range in ionization parameter. 

We fit the SISO and SIS1 data simultaneously and obtain a statistically acceptable fit with 
Xr = 1.10 with 290 degrees of freedom. The data in comparison to the model are plotted in 
Figure 3. In addition to the continuum parameters, this model contains normalizations of each ion 
line flux, which are proportional to their ion line emission measures. Figure 4 shows the derived 
empirical DEM distribution assuming an emission width of Alog£ ~ 0.5 for helium- like ions and 
Alog^ ~ 1 for the hydrogenic species for two sets of assumed chemical abundances; (1) solar 
photospheric abundances (Anders & Grevesse 1989) and (2) abundances derived by Oliva, Marconi, 
&; Moorwood (1999) for a knot located approximately 15" from the nucleus in the direction of the 
[O III] cone. We adopt solar abundances for elements that are not listed in Oliva, Marconi, & 
Moorwood (1999). Given the assumptions that we have made in the construction of the empirical 
DEM distribution, we estimate the points drawn on Figure 4 to be uncertain by up to a factor of 2 
3. 

With the exception of the DEM points derived from O VII and O VIII at log£ = 1.4 and 
log£ = 1.9, respectively, the shape of the empirical DEM is nearly linear with evidence of a slight 
decrease in DEM with increasing log£. The oxygen ions have bright features below E ~ 1 keV 
and the derived emission measures are extremely sensitive to the foreground column density. An 
increase in the column density by 50% can be offset by an increase in the emission measure of 
O VII by a factor of ~ 10 and O VIII by a factor of ~ 6. Keeping this in mind, the derived DEM 
distribution, to first order, is nearly linear. 

We then adopt a reverse fitting procedure where we parametrize the model spectrum of the 
highly ionized component according to powerlaw DEMs of the form, d(EM)/dlog£ oc £ 7 , with 
minimum and maximum cutoff ionization parameters, £ m i n and £ maX ) respectively. This procedure 
bypasses the ambiguities in the assumed temperatures of formation and correctly takes account 
of the temperature gradient in each ionization zone. In addition to the K-shell ions that we had 
initially used in the empirical fit, this time we also include line and RRC emission from the iron L 
species, which are expected to contribute ~ 10% to the total recombination emission flux. There are 
four free parameters in this component, namely, the slope 7 of the DEM as a function of ionization 
parameter, the normalization, which is proportional to the total emission measure EM, low-£ 
cutoff £ m i n , and high-£ cutoff £ max - Unlike the CIE case, we are able to centroid the fluorescent line 
energies of all the elements that we list in Table 2, and we find a statistically acceptable fit with 
Xr = 1.12 for 299 degrees of freedom. The data in comparison to the model is shown in Figure 5 
and the best-fit parameters are listed in Table 3. We also derive confidence contours for two pairs 
of interesting parameters; (£ m inj Cmax) and (7, EM), as shown in Figures 6 and 7. In summary, 
we derive 90% confidence ranges of 7 = — 0.07!q'i4 an d EM = (5.5 ± 0.7) x 10 63 cm~ 3 for two 
interesting parameters, and log£ m i n = 0.9ll^, 8 and log£ max = 3.53l '2Q, again, for two interesting 
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parameters. The best-fit model spectrum with the various components is shown in Figure 8. 

The lower limit on £ m i n cannot be well-constrained simply because the X-ray spectrum of the 
Circinus Galaxy is absorbed through a moderately high column density (Njj ~8x 10 21 cm -2 ), and 
also because the ASCA SIS detectors are not sensitive to X-rays below ~ 0.5 keV, where discrete 
recombination emission from regions with log^ ~ 1 dominates, for example, from carbon and 
nitrogen. The upper limit on £ max is constrained mainly by the Lya flux of Fe XXVI. 

4. Discussion and Implications 

We have shown that the soft X-ray emission line spectrum of the Circinus Galaxy is compatible 
with both coronal and photoionization equilibrium, and the excitation mechanisms cannot be spec- 
troscopically distinguished without any ambiguity. In this section, we extract physical parameters 
implied by each model and discuss their validity in relation to observations at other wavelength 
bands. 

According to recent ROSAT and ASCA observations, nearby star-forming galaxies show X- 
ray spectra that can be well described by a combination of a soft optically thin thermal plasma 
with a temperature of kT ~ 0.6 keV and a hard continuum component with either a photon 
index of V ~ 1.7 or a bremsstrahlung temperature of kT ~ 5 keV (e.g., Ptak et al. 1999, and 
references therein). The collisional ionization interpretation for the Circinus Galaxy gives a best-fit 
temperature of ~ 0.69 keV and a luminosity of few x 10 40 erg s _1 in this component, which are 
typical values for starburst galaxies. Note that the temperature of the plasma is well-constrained 
and that the data do not allow the presence of gas with temperature lower than kT ~ 0.60 keV at 
the 99% confidence level. It is interesting to compare this with what Oliva et al. (1994) have found 
using visible and near infrared observations of the coronal emission lines in the Circinus Galaxy. In 
order to reproduce the observed spectrum in the context of a powerlaw DEM plasma in CIE, they 
require a sharp cutoff in the DEM at a temperature of kT max ~ 0.1 keV. This would imply that 
there is a gap in the emission measure distribution in the range 0.1 < kT < 0.6 keV, which has no 
obvious physical interpretation. 

In PIE, the assumption of plasma heating by a point source of X-radiation places an additional 
constraint on the plasma structure (e.g., the heating rate as a function of £), which can be then used 
to derive the geometrical distribution of the gas. In the PIE interpretation of the spectrum, have 
shown that a nearly flat DEM distribution is consistent with the data. One simple way of producing 
a flat DEM curve is to adopt a radial density profile of the form, n e (r) oc r -3 / 2 . This is perhaps 
the simplest configuration but, by no means, unique. More generally, a powerlaw density profile 
produces a powerlaw DEM curve. It is interesting to note, however, that a conical homogeneous 
medium with a powerlaw density profile as a function of radius always produces a powerlaw DEM 
distribution; i.e., a simple calculation shows that for a density profile of the form, 

n e (r) = n e0 r^f , (1) 
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where we define n e o to be the density in cm 3 at a distance of 1 pc from the central source, the 
DEM distribution can be written as, d(EM) = AO, n 2 . r 2 dr, or 

d ( EM ^ - l n 10 A ^ \ n n (uL )3-2/3l 1/(2-/3) , (3-2/3)/(/3-2) (2) 

where /x is the ratio of the electron to proton number densities (fi = 1.2, for an X-ray emitting 
plasma) and AO, is the opening solid angle of the photoionized region. From this equation, we can 
identify the slope of the DEM curve 7 to be (3 — 2(3)/ {(3 — 2). Our inferred range in 7 corresponds 
to (3 = 1.48^q'q4. One can also see that (3 = 3/2 corresponds to a flat DEM curve. 

The fact that (3 is relatively well-constrained is due to the fact that the slope of the DEM curve 
is extremely sensitive to the value of (3. As mentioned earlier, a (3 = 3/2 density profile produces a 
flat DEM curve, which produces emission lines and RRC from a wide range of ionization stages. A 
(3 = 2 profile, on the other hand, corresponds to a single-^ zone, since both the ionizing flux and 
the density decay like oc r~ 2 . In this case, only a small number of charge states can exist, and the 
resulting spectrum is more sparse. 

Our derived values for the parameters shown in Table 4 depend on the assumed metal abun- 
dances. If we assume that the abundances are Az e relative to solar, the DEM curve will be 
uniformly shifted by a factor of Ag 1 and the inferred total emission measure will be EM = 
5.5 x 10 63 Ag 1 cm~ 3 . The parameters, (3, £mim an d £max are unaffected by a uniform shift in 
abundances, since they characterize only the shape of the DEM curve and not its normalization. 
The exact dependences of these parameters on the abundance of an individual element, however, is 
more complex. For example, an underabundance of oxygen and neon, as found by Oliva, Marconi, 
& Moorwood (1999), may result in a steeper DEM slope {(3 > 1.5) or a decrease in £ m j n in order to 
increase the emission measure at low-£ and compensate for the observed flux at low energies. 

Although the density parametrization of equation 1 is simplistic, chosen for the photoionization 
interpretation of the spectrum, there is empirical evidence that narrow-line regions in many Seyfert 
galaxies, in fact, do exhibit this type of behavior. For (3 = 2, the density scales like r~ 2 , and, 
hence, the ionization parameter is constant with radial distance. This has been observed in some 
systems that show correlations between widths of optical emission lines and their critical densities 
for collisional de-excitation. Spectral transitions with higher critical densities generally exhibit 
broader line profiles. There are also cases where the observed widths are higher for transitions 
in ions with higher ionization potential, which corresponds to a scenario with (3 < 2 (e.g., Pelat, 
Alloin, & Fosbury 1981; Filippenko & Halpern 1994; De Robertis & Osterbrock 1984). According 
to De Robertis & Osterbrock (1986), values of (3 lie between and 2, while more Seyfert 2 galaxies 
have (3 closer to 2. A (3 > 2 behavior (i.e., £ increasing with distance) is not observed in any Seyfert 
galaxy. 

Using the empirical density profile, we can estimate several physical parameters from the 
results of the DEM fit. From the best-fit values for (3 and the total EM, which is proportional 
to the prefactor in equation 2, Ar2(n e o-i^~ 2/3 ) 1 /( 2_/3 - ) , we estimate the density profile of the X-ray 
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emission region to be, 

n e (r) = 4.0 x 10 2 AJT ' 52 r" 1 ' 48 (cm" 3 ), (3) 

where L x ^ 2 ) is the X-ray luminosity in multiples of 10 42 erg s _1 . We also saw that the upper 
and lower limits on the ionization parameter were log£ 

max ~ 3.5 and log£ m ; n ^ 0.9, respec- 
tively. From the upper limit, we can estimate the distance from the central nucleus to the inner- 
most X-ray emission line region that is directly visible from the observer. We estimate the inner 
distance to be r in ~ 0.011 AS1 L 2 x ^ pc. Similarly, the lower limit on the outer distance is 
?out ~ 1200 AVt £ 2 ( 42 ) pc We can see from the [O III] image that the opening half-angle of the 
one-sided ionization cone is ~ 30°, which corresponds to a solid angle of AQ ~ 0.8. Assuming 
that the X-ray emission lines are produced in the inner region close to the base of this ionization 
cone, we estimate its size to be ~ 900 pc or approximately 50" in the sky. This is qualitatively 
consistent with the ROSAT HRI image that shows a slight spatial extent (< 30") approximately in 
the direction of the [O III] cone. 

As a cross-check on the derived parameters, we calculate the electron scattering optical depth 
through the emission line medium and compare it with the observed scattered continuum flux. As 
listed in Table 3, the unabsorbed luminosity of the scattered continuum is L s x cat ~ 8.4 x 10 39 erg s -1 , 
which corresponds to a Thomson optical depth of r es = 0.084 (AQ £ x (42)) -1 - Integrating the density 
profile (equation 3) from r m out to r out yields an optical depth of, r cs = 0.014 (AilL^^) -1 . We 
note that regions beyond r out do not contribute significantly to the optical depth. This implies 
the presence of additional scattering regions in the circumsource medium, which must be either 
more highly ionized or less ionized than the X-ray emission line region. The additional column 
density required is Nh (Af2/47r) ~ 1 x 10 23 cm~ 2 , which may be produced in the region that 
absorbs the transmission component observed in the hard X-rays. The measured column density 
of this component is Nh ~ 4 x 10 24 cm~ 2 (Guainazzi et al. 1999), and thus a covering fraction of 
(A$7/47r) ~ 0.03 can produce most of the observed soft continuum flux. An optical depth of r es ~ 
0.01 through the photoionized gas amounts to a total column density of Nh ~ few x 10 21 cm~ 2 , 
which justifies our assumption that the primary radiation is optically thin to continuum absorption. 

It is interesting to note that a density profile of the form n(r) ~ r~ 3 / 2 is expected for Bondi 
accretion in a spherically symmetric gravitational potential of a point source. For a central black 
hole mass of M§ (in multiples of 10 6 M ), the gas is nearly in free-fall within an accretion radius 
of r acc ~ 100 Mq ^W4) P c > where r oo ( 4 ) is the ambient gas temperature (in multiples of 10 4 K) 
far away from the central compact object. Using the density profile of equation 3 and a free-fall 
velocity profile (i.e., v oc r -1 / 2 ), the mass inflow rate can be expressed as, 

M inflow = 8 x 10- 4 Aft ' 48 L-° 42 04 Ml 12 (M yr" 1 ). (4) 

The mass accretion rate implied by the observed X-ray luminosity is, 



Maccretion = 2 x 10 4 L x(42) rj Q \ (M© yr *), 



(5) 
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where 770.1 is the energy conversion efficiency in multiples of 0.1. Since the accretion radius is on the 
order of or larger than the size of the X-ray emission line region, the outer region may not necessarily 
be inflowing at the free-fall velocity, especially since photoelectric heating by the AGN continuum 
radiation increases the local sound speed of the surrounding gas. Nevertheless, these mass flow rates 
are roughly consistent with each other and we suggest that the inner part of the X-ray emission 
region is fueling the central engine. Projected velocity shifts of up to ~ 230 km s _1 have also been 
detected near the base of the [O III] and Ha ionization cones, which may be interpreted as either an 
outflow or inflow relative to the nucleus (Veilleux & Bland-Hawthorn 1997; Elmouttie et al. 1998). 

We suggest that the region that produces the observed X-ray recombination lines is the con- 
fining medium of the so-called coronal line regions (CLR; Oke & Sargent 1968). Using the surface 
brightness of [Fe XI], (Oliva et al. 1994) finds the CLR to have a size ~ 0.5" (diameter of ~ 10 pc) 
and a gas density of n e ~ 250 f^ 1 ^ 2 cm" 3 , where / is the volume filling factor. The temperature 
is assumed to be kT ~ 3 x 10 4 K, a typical value for a photoionized CLR. Using the density profile 
inferred from our fit, we estimate the gas density at a distance of ~ 10 pc to be n e ~ 10 cm~ 3 . 
According to our model, the ionization parameter at this distance is log£ ~ 2.5 with a temperature 
of kT ~2x 10 6 K, and thus the pressure is n e T ~ 2 x 10 7 . In order for the two regions to be in 
pressure equilibrium, the filling factor of the CLR must be / ~ 0.1. This implies the density of the 
[Fe XI] region to be n e ~ 800 cm~ 3 , which is still much lower than the critical density of the upper 
level and, therefore, the line will not be diminished. 

The model that we have suggested can be further tested with high spectral and spatial res- 
olution grating observations with the Chandra observatory. Our model predicts that ~ 70% of 
the total flux will be spatially unresolved, containing the reflection component, as well as a large 
fraction of the scattered soft continuum component. The remaining ~ 30%, which is dominated 
by recombination line emission, will be extended. Specifically, recombination lines from \ow-Z ions 
(e.g., O, Ne, and Mg) will be extended by ~ 10" and those from higher- Z ions (e.g., Si, S, Ar, 
Ca, and Fe) will have a spatial extent of < 5". In addition to the spatial distibution of the X-ray 
emission line region, measurements of velocity shifts will allow us to constrain the kinematics and 
determine whether the confining medium is flowing towards or away from the central nucleus. 

This work has benefited from useful discussions with Andy Rasmussen concerning the ASCA 
SIS calibrations. We have made use of data obtained through the High Energy Astrophysics Science 
Archive Research Center Online Service, provided by the NASA/Goddard Space Flight Center. We 
thank the ASCA team and the HEASARC Online Service for their support. This work was supported 
under a NASA Long Term Space Astrophysics Program grant (NAG5-3541). D. A. L. was supported 
in part by a NASA Long Term Space Astrophysics Program grant (S-92654-F). Work at LLNL was 
performed under the auspices of the U. S. Department of Energy, Contract No. W-7405-Eng-48. 
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Fig. 1. — The ASCA SIS data and the best-fit CIE model. The fit is statistically acceptable with 



XT 



1.21 for 305 degrees of freedom. Only SIS0 data is shown for clarity. 
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Fig. 2. — Two parameter confidence ranges (kT vs. total EM) for the highly ionized component 
in the CIE fit. Solar abundances are assumed. 
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Fig. 3. — The best-fit recombination cascade model compared with the SIS data (Xr = 1.10 with 
290 degrees of freedom). Again, only SISO data is shown for clarity. 
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Fig. 4. — An empirical DEM distribution derived from K-shell oxygen, neon, magnesium, silicon, 
argon, and iron for two sets of assumed chemical abundances; solar photospheric abunandes of 
Anders & Grevesse (1989) and abundances derived by Oliva, Marconi, & Moorwood (1999) for a 
knot located near the nucleus. 
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Fig. 5. — The best-fit powerlaw DEM PIE model in comparison to the data. Again, the fit is 
acceptable with Xr = 1-12 for 299 degrees of freedom. 
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Fig. 6. — Two parameter confidence ranges for the slope of the DEM curve, 7, versus the total 
emission measure, EM, for the PIE DEM fit. 
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Fig. 7. — Two parameter confidence ranges for the lower and upper cutoff ionization parameters 
(log£min vs. log«e max ) for the PIE DEM fit. 




Fig. 8.~ 



The best-fit PIE DEM model spectrum for the parameters shown in Table 4. 
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Table 1. Summary of Fluorescent Line Intensities 



CIE Interpretation PIE Interpretation 51 



Element Line Energy b (keV) Line Intensity Ax 2 Line Intensity Ax 2 Reflection 01 



Si l-74to3o °- 08 L2 °- 23 18 - 6 °- 18 

S 2.34+°;°^ 0.26 17 0.30 22.4 0.17 

Ar 3.05±g;gf 0.062 3.4 0.08 9.2 0.066 

Ca 3.73lg;5o °- 018 °- 06 °- 067 L3 °- 051 

Fe (Kq) 6.49±g;gi 3.1 520 2.9 540 3.0 

Fe (K/3) 7-18lao5 0.42 59 0.40 49 0.38 

Ni 7.5612-J2 0.038 1.3 0.013 3.5 0.10 



a Derived fluorescent line intensities when the highly ionized line emission component is 
modeled as plasmas in CIE and PIE. 

b Measured line centroids in the PIE fit, which are fixed in the CIE fit. 
°In multiples of 10~ 4 photons cm~ 2 s _1 

d Line intensities expected from pure cold reflection normalized to the average observed 
Fe Ka flux. 
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Tablc 2. Summary of the Collisional Ionization Equilibrium Fit 



Parameter Value 

N H 1.4 x 10 22 cm" 2 

r a 1.5 

L soft b 3 x 10 40 erg s" 1 

fcT keV 

EM (1.12lg;^) x 10 63 cm" 3 



a Photon index of the scattered and 
reflected continuum components, which 
are set equal in the fit. 

b Derived luminosity in the soft con- 
tinuum component. 
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Table 3. Continuum Parameters for the Photoionization Equilibrium Fit 



Values 



Parameter Empirical DEM 

N H (10 22 cm" 2 ) 0.87 0.89 

r a 1.8 1.6 

L soft (erg s" 1 ) 1.2 x 10 40 8.4 x 10 



a Photon index of the scattered and re- 
flected continuum components, which arc 
set equal in the fit. 
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Table 4. Summary of the Photoionization Ionization Equilibrium Fit 



Parameter Value a 
P 1.48±^ 

£ q cq+0.32 

?max o.JO_ q 20 

EM 5.5 ± 0.7 x 10 63 cm" 3 

a Errors correspond to 90% confi- 
dence ranges for two interesting pa- 
rameters as described in the text. 



